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Abstract 

This article is intended for physics educators and students at school 
and undergraduate level. It is used at our department to introduce stu- 
dents to simulation and offer a guide in using statistics in physics. The 
simulation code was created using Matlab, and was given a friendly 
interface with a Labview module. 



1 Introduction: 

During our teaching for undergraduates of physics we saw a need for an 
experiment that can demonstrate the close ties between experiments and 
statistical analysis. We also felt the need of introducing the idea of simulating 
a physics experiment for the huge and increasing role of simulation in the 
different fields of science. We used to have an experiment targeting the 
Poisson distribution and Radioactivity, the experiment was performed with 
a Geiger counter, so it was a good choice for combining the two ideas. The 
simulation code was created with a Matlab script and was then given a 
friendly interface with the help of a Labview module. For running the code 
you need to have Matlab installed on your system. 
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2 Theory: 



The simulation depends on the Normal number generator available through 
the Matlab. To simulate the decay behavior of a radioactive sample having 
N a nuclide, and decay constant A during a time At, N Q uniform random 
numbers were generated, with each number representing a potential decaying 
nuclide. To demonstrate this note that the probability of having a random 
number x between a and b such that x Q > x > b is given by 



P(x > x ) = f 

J X 



T^—dx; 
b — a 

(b - x ) 



(1) 



b — a 

setting a = and b = 1, in our case we get 

P(x > x ) = 1 - x . (2) 

Noting that the probability of a single nuclei to decay in a time At is XAt, 
or 



A At = 1 — x , hence 
x = 1 - XAt. (3) 

i.e any uniform random number grater than x Q represent a decaying nuclide. 
The Matlab uniform random number generator is used to create iV random 
number and each is checked with condition equation© to identify it as a 
non-decaying or decaying nuclei which in such a case counted. 

The time At is controlled trough the Labview interface. In the simulation 
you are given a number of different samples to choose among each represent 
a different decay constant. You also have three choices of N such that for 
small N a = 10 6 , medium N Q = 5 x 10 6 and large N Q = 50 x 10 6 . The above 
choices were made for the limited computational power of the system used. 
The number of trials to be simulated can be controlled and up to a thousand 
trials can be simulated. The result of each trial is recorded at a .txt file of 
your choice. 
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3 Application: 



The Labview is started, and after making your choices of the different vari- 
ables press the start bottom. The result of each trial will be recorded at a 
file, which can then be read and the data treated with a data analysis tool. 
An example of the simulation is the histogram below created from 851 trials 
for sample E with medium number of nuclei for a time of 60min for each 
trial. Using all this information the half life of the sample can be calculated 

N — N (l — e- XAt ) . (4) 
Using the value of N collected from the simulation, A can be found as 

79.189 = 5 x 10 6 (l - e~5«W) , 

Ti = A.990years. 
Compare this to the value set in the code. 



Note: 

If you would like to have the simulation please send me an email and we will 
be happy to send it. 
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Sample E 




Counts 
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4 The Simulation: 



In this section we include the code and a snapshot of the Labview interface. 
Matlab Script 

1: format long g 2: if n==0; 

3: x= l-(log(2)/(0.5*24*3600))*l*60*t; 

4: elseif n==l; 

5: x= l-(log(2)/(10*24*3600))*60*t; 
6: elseif n==2; 

7: x= l-(log(2)/(160*24*3600))*60*t; 
8: elseif n==3; 

9: x=l-(log(2)/(2*365*24*3600))*60*t; 
10: elseif n==4; 

11: x=l-(log(2)/(5*365*24*3600))*60*t; 

12: end 

13: if c==0; 

14:d=1000000; 

15: elseif c==l; 

16: d=5000000; 

17: elseif c==2; 

18:d=50000000; 

19: end 

20: z=0; 

21: for i=l : d 

22: y=rand; 

23: if yix 

24: z=z+l; 

25: else 

26: z=z; 

27: end 

28:end 

29: z 
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